Water transport in protoplanetary disks and the hydrogen isotopic composition of 

chondrites 



Emmanuel Jacquet a ' b '*, Francois Robert b 



" Canadian Institute for Theoretical Astrophysics, University of Toronto, 60 St Georges Street, Toronto, ON M5S 3H8, Canada 
b Laboratoire de Mineralogie et de Cosmochimie du Museum, CNRS & Museum National d'Histoire Naturelle, UMR 7202, 57 rue Cuvier, 75005 Paris, France. 



Abstract 



m 
o 

(N 
C 

Hi 

m 

(N 



PL, 

m 

Oh 
i 

o 

» 

03 



The D/H ratios of carbonaceous chondrites, believed to reflect the hydrogen isotopic composition of water in the inner early solar 
system, are intermediate between the protosolar value and that of most comets. The isotopic composition of cometary water has been 
accounted for by several models where the isotopic composition of water vapor evolved by isotopic exchange with hydrogen gas 
in the protoplanetary disk. However, the position and the large range of variation of the distribution of D/H ratios in carbonaceous 
chondrites have yet to be explained. In this paper, we assume that the D/H composition of cometary ice was achieved in the disk 
building phase and model the further isotopic evolution of water in the inner disk in the classical T Tauri stage. Reaction kinetics 
compel isotopic exchange between water and hydrogen gas to stop at ~500 K, well inside the snow line. However, the equilibrated 
water can be transported to the snow line (and beyond) via turbulent diffusion and consequently mix with isotopically comet-like 
water. Thus the competition between outward diffusion and net inward advection established an isotopic gradient, which is at the 
origin of the large isotopic variations in the carbonaceous chondrites and other water-bearing objects accreted in the protoplanetary 
disk. 

Under certain simplifying assumptions, we calculate analytically the probability distribution function of the D/H ratio of ice 
accreted in planetesimals and compare it with observational data. The distribution is found to essentially depend on two parame- 
ters: the radial Schmidt number Scr, which ratios the efficiencies of angular momentum transport and turbulent diffusion, and the 
range of heliocentric distances over which currently sampled chondrite parent bodies were accreted. The minimum D/H ratio of 
the distribution corresponds to the composition of water condensed at the snow line, which is a function of both the composition 
of equilibrated water having diffused from hotter disk regions and the efficiency of this outward transport as measured by Sc«. 
Observations constrain the latter to low values (0.1-0.3), which suggests that turbulence in the planet-forming region was hydrody- 
namical in nature, as would be expected in a dead zone. Such efficient outward diffusion would also account for the presence of 
high-temperature minerals in comets. 
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1. Introduction 

The isotopic composition of hydrogen (as expressed by the 
D/H ratio) is a valuable tracer of the origin of water in the so- 
lar system (Robert 2006). Primitive meteorites (chondrites), 



presumably the building blocks of terrestrial planets, allow a 
glimpse at the D/H ratio of water in the protoplanetary disk 
that surrounded our Sun 4.57 Ga ago. Indeed, many chondrites, 
in particular carbonaceous chondrites, contain clays formed 
through alteration of originally anhydrous silicates by water 
presumably incorporated as ice along with rock during accre- 
tion ( |Brearley|20"0"3"t|Ghosh et aZ.|2006| l. 

The measured D/H ratios of bulk carbonaceous chondrite^] 
(see Fig. [5]), g enerally thought to reflect that of accreted water 
(but see Alexander et al. 1 |2012b| , span a range of 120 x 10 



to 230 x 10 6 (excluding CR chondrites). The distribution, 
which is skewed to heavy isotopic compositions, has a mean 
(156 + 3) x 10" 6 close to the (149 ± 3) x 10~ 6 estimated for the 
bulk Earth ( Lecuyer etal.\\ 998 1 — consistent with a chondritic 
source for terrestrial water. The D/H ratios of carbonaceous 
chondrites is systematically lower tha n that exhibited by m any 



Oort-cloud comets ((296 ± 25) x 10"; Hartoghef al. 2011 1, al- 



though D/H ratios of (161 + 24) x 10 6 and (206 + 22) x 10" 6 
have been determined for Jupiter-family comet Hartley 2 and 
Oort-cloud comet Garradd, respectively (Hartogh e t q/.||20lT 



Bockelee-Mor van et a/.|2012) . The composition of interplane- 



tary dust particles, though broadly similar to that of carbona- 
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ceous chondrites (e.g. Engr and and Maurette| 1998 Bradley 
2005), has a significant tail extending to cometary values and 
beyond. Both chondritic and cometary domains of variation 
are markedly distin ct from both th e estimated protosolar value 
Geiss and Gloeck"Ier||2003 1 



((20 ± 3.5) x 10" 



dominated by 



the isotopic composition of hydrogen gas, and the high values 
(D/H > 10 3 ) determined in molecular clouds for molecules 
other than H2, consistent with predictions from ion-molecule 
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Figure 1: Cartoon of the model investigated here. We distinguish between 
(comet-like) "cometary water" (cyan) and "equilibrated water" (pink) that has 
undergone isotopic exchange with the hydrogen gas inside the "equilibration 
line" (R le ac)- Water is gaseous inside the snow line (S con( j) and solid beyond 
(as symbolized by the cloud and snowflakes symbols, respectively). Accretion 
of water as ice occurs beyond the snow line. Arrows symbolize motions of the 
gas. 



reactions (see e.g. Robert 2006 1. However, high D/H values 
> 2 x 10 3 have been reported at the micrometer scale in clays 
of the Semarkona chondrite (Pian i et a/.|2012||Piani|2012l ) and 
may reflect the composition of pristine interstellar ice grains. 

In the protoplanetary disk, the D/H ratio of water evolves 
mainly through isotopic exchange with the hydrogen gas 
(whose composition remains essentially fixed at the protosolar 
value because it contains the bulk of the H of the system), that 
is, via the reaction: 



HDO + H 2 H 2 + HD 



(1) 



The equilibrium fractionation factor (D/H) H ,o/(D/H) H2 is 
unity at high temperatures (> 1000 K), but increases with de- 
creasing temperature ( Richet et al.\\911 \. 



Using one-dimensional disk models, Drouart et al. ( 1999 1 
showed that reaction ([T} was unable to account for the chon- 
dritic or cometary D/H ratios if water was assumed to have 
formed in the protoplanetary disk with the protosolar D/H 
value. Indeed reaction kinetics would have been prohibitively 
slow in those cold regions of the disk where sufficient D en- 
richment would have been predicted by thermodynamic equi- 
librium. On the other hand, |Drouart et al.\ ( | 1999) 1 found that 
cometary compositions could be obtained if water was inher- 
ited from the protosolar cloud (with heavy D/H ~ 10~ 3 ). Water 
would have evolved toward isotopically lighter compositions 
in the inner disk and diffused in the colder regions owing to 



turbulence, a result also obtained by Mousis et al. ({2000 ) and 
Hersant et al. (|2001|l. Transport and mixing would have been 



particularly efficient in the early stages of disk evolution, if the 
disk was first built compact and then expanded because of tur- 
bulence (Yang et al. 2012). Such a picture is also consistent 



with the presence of crystalline silicates in comets (Bockelee- 
|Morvanef q/.|2002| l. 

While these previous works focused on reproducing D/H 



value observed in comets, implications on the composition of 
water in the inner solar system, and in particular the fairly large 
domain of variations of the D/H ratio in chondrites have yet to 
be worked out in this framework. This is the purpose of this 
paper. 

In this paper, we consider a disk model pertaining to the clas- 
sical T Tauri phase, that is, after infall has ceased and the disk 
radius has become large compared to the heliocentric distances 
of interest. This is the epoch where chondrite accretion is be- 
lieved to have taken place. At that time, the D/H ratio of water 
ice in the outer disk is assumed to have been already set and 
homogenized at the value of ~ 300 x 10~ 6 measured in most 
comets as a result of early disk processes (this would corre- 

( |20T2] l, 



spond to the plateau in the simulation of Yang et al. 



whose level varies little after infall has stopped). Close to the 
Sun, isotopic exchange between (gaseous) water and hydrogen 
gas (reaction (|T]i) is still going on. However, beyond a certain 
heliocentric distance (corresponding to a temperature of ~500 
K, as we shall argue later), the kinetics of this reaction are so 
slow that essentially no isotopic exchange occurs there. Due 
to gas turbulence, however, some of the water equilibrated with 
H 2 sunward of this heliocentric distance will be transported out- 
ward and will reach the snow line, where ice condenses. Water 
beyond the snow line — sampled by chondrites — will thus be a 
mixture of two components: equilibrated water having diffused 
from the hot inner regions of the disk, and cometary water drift- 
ing inward from the outer disk. 

A radial gradient in D/H ratio is thus established by the com- 
petition between outward diffusion and inward advection which 
sets the proportions of equilibrated water and cometary water as 
a function of heliocentric distance. The range of D/H ratios of 
water-bearing bodies thus extends from cometary compositions 
down to the composition of water condensed at the snow line. 
The composition sampled at the snow line is hence intermediate 
between cometary values and those expected from equilibration 
with H2. This schematic description provides a framework to 
account for the fairly large variations of the D/H ratio in water 
of solar system bodies. A cartoon of the scenario is presented 
in Fig. [1] 

The goal of this paper is to calculate analytically the distri- 
bution of D/H ratios of chondritic water that results from this 
picture and compare it to the observational data, in order to con- 
strain the model parameters. In particular, the calculation will 
be constrained by the asymmetrical shape of the distribution 
and its position relative to cometary values. We stress that here 
the cometary value is a starting, observationally given param- 
eter of the model. Our model does not aim at reproducing the 
composition of comets in contrast to the studies of Drouart et al. 



(1999); Mousis et a l.\ ( 12000) 1; [Hersant et aZ. | ( [200 l| l; | Yang et al. 
I20T2I ; rather, in complementarity to those, it focuses on chon- 
drite parent bodies, that is, the inner regions of the solar system. 
Also in complementarity to these numerical studies, our work, 
in being analytic in nature, allows us to nail down, under cer- 
tain simplifying assumptions, the relevant parameters (essen- 
tially two) that govern the distribution of chondritic D/H ratios, 
namely: the range of heliocentric distances sampled by chon- 
drites and the radial Schmidt number — which essentially ratios 



2 



the efficiencies of angular momentum transport and turbulent 
diffusion. We will find that low values of the radial Schmidt 
number yield D/H variations consistent with observations. The 
paper is organized as follows: We describe the model assump- 
tions in Section|2] present and discuss the results in Sections [3] 
and|4] respectively. In Section [5] we conclude. For the sake of 
clarity, specific derivations are deferred to appendices. 

2. Modeling 

In this section, we introduce our notations and modeling prin- 



ciples. We successively consider the disk model (Section 2.1 



the transport of water (Section 2.2 1, its D/H ratio (Section 2.3 1 



and our prescription for accretion and delivery to Earth (Section 
2.1. The disk 

We consider an axisymmetric, vertically isothermal turbulent 
disk. The disk is assumed to have a stationary surface density 
profile in its inner regions such that its (radially constant) mass 
accretion rate M is given by: 



M = -2nRLu R = ?>7tY.a 



(2) 



with R the heliocentric distance, 2 the surface density, ur the 
vertically averaged radial velocity of the gas, a the vertically 
averaged turbulence parameter, £2 the Keplerian angular ve- 
locity and c s - y/ksT/m the isothermal sound speed — where 
kg and m are the Boltzmann constant and the mean molecular 
mass, respectively, and T the temperature. The steady-state ap- 
proximation is expected to hold so long the viscous evolution 
timescale, 



tms(R) = 



R 1 



3ac 2 s /n 



= 0.05 Ma R 



/300K\/10 



(3) 



is shorter than the disk's age in the regions of interest, 
ilarly to |Drouart et al. \\999) and IHersant et al. 



also IHartmann et al. 



the mass accretion rate to evolve in time a; 



Sim- 
( pOOT) — see 



998]i~and [Chambers! ( |2009| >— , we take 



M = 



M 



(1 + t/fo) 



3/2 



(4) 
= 



where Mo is the mass accretion rate at some initial time t 
and to an evolution timescale (to > f v ; s (/?)). 

We adopt the same temperature prescription as in Jacquet 

\eTaT\$2W2\ : 



max 



1287T 2 cr SB fe B Q; 



M 2 Q 3 



1/5 



JtToR^ 2 



(5) 



with k the specific opacity, cr SB the Stefan-Boltzmann constant, 
Tq = 280 K, fx a dimensionless constant parameter (for which 



2 This is the evolution expected fro m the self-sim ilar solution of the conti- 
nuity equation for ac'JCl oc R (see e.g. |Garaud|2007| as is the case here in the 



we will adopt a fiducial value of 0.5) and R AV = R/(l AU). 
This prescription essentially states that inner disk regions are 
dominated by the dissipation of turbulence ("viscous heating"), 
which decreases with decreasing accretion rate, whereas the 
outer disk regions are dominated by reprocessing of solar ra- 
diation. Hence, the snow line (at heliocentric distance ^ con d), 
which corresponds to a fixed temperature of T COB d=170 K, re- 
cedes toward the Sun with the passage of time. 

For simplicity, we will assume that a and k are constant 
throughout the disk's extent and evolution. 

2.2. Water transport 

We now turn to the transport of water. We distinguish be- 
tween nebular and accreted water. Nebular water is water dy- 
namically coupled to the gas, whether as water vapor or ice- 
bearing grains, and behaves as part of the disk. It is predom- 
inantly gaseous inside the snow line and solid outside it. Ac- 
creted water refers to water retrieved from the gas by incorpo- 
ration in planetesimals or comets. 

We first focus on the dynamics of nebular water. Since it be- 
haves as a passive contaminant in the gas, the evolution equa- 



tion of its surface density 2h,o reads (e.g. Ciesla and Cuzzi 
|2006) : 



d£ H2 o }_d_ 
dt + RdR 



6r ^3R 



£h 2 o 

= S coll 



(6) 



with S co n the sink term resulting from collisions — both co- 
agulation and shattering processes (contributing positively and 
negatively to it, respectively) — , and Sr a dimensionless num- 
ber (of order a) parameterizing turbulent diffusion. We neglect 
any gas-solid drift induced by gas drag (we shall return to this 



assumption in Section 4.2 1 



If coagulation and shattering can be ignored, it is straight- 
forward to see that 2h,o x £ is a solution of the equation. 
This should remain a reasonable approximation for coagula- 
tion timescales longer than ? V i S . There is evidence suggesting 
that accretion was indeed quite inefficient/slow in the early so- 
lar system: age dating of chondrite components and studies of 
the thermal evolution of their parent bodies are consistent with 
accretion being a protracted process on a few-Ma timescale 



(e.g. Villeneuve et al. 2009 Grimm and McSween 1989 Kleine 



\et a/.|2008||Connelly et aZ.|20122 and the total mass of the cur 
rent planetary system (see e.g. Hayashi 1981) is one order of 
magnitude lower than the nonvolatile content of disks at the end 
of the infall phase (e.g. Yan g and Ciesla|2012 1 assuming a solar 
metallicity. We shall thus assume that accretion was inefficient 
and ignore its feedback on the dynamics of nebular water so 
that ee,o = £h,o/£ can ^ e considered constant (we shall return 
to this issue in Section[4]i. 

2.3. Hydrogen isotopic composition of water 

In this paper, we assume that the D/H ratio of water is 
dictated by isotopic exchange between originally isotopically 
heavy (comet-like) water and isotopically light hydrogen gas 



outer disk. 



(e.g. Drouart et al. 1999 Mousis et al. 2000 Hersant et al. 



3 



2001 1. To make the problem analytically tractable, we circum- 



vent an accurate treatment of reaction kinetics by schematically 
distinguishing between two types of water: cometary water is 
water that has never experienced temperatures in excess of a 
"reaction temperature" r reac during the classical T Tauri stage 
of the disk, and is assumed to retain a relatively heavy D/H ratio 
denoted (D/H)/,; equilibrated water, on the other hand, is water 
that did experience temperatures above r reac and therefore has a 
D/H ratio set to a fixed value (D/H); through isotopic exchange 
with hydrogen gas (a justification for this treatment based on 
reaction kinetics is provided in Appendix B I. T relLC , (D/H)/ and 
(D/H)/, are fixed parameters of the model and are assigned val- 
ues of 500 K, 40 x 10~ 6 and 300 x 10~ 6 for them, respectively 
(see Appendix B i. The value of 40 x 10~ 6 corresponds to the 



maximum D/H ratio that can be reached by water in equilib- 
rium with molecular hydrogen. Indeed, although enrichment in 
D could be in principle higher at lower equilibration temper- 
ature, slow kinetics of the isotopic exchange reaction prevent 
this equilibrium from being actually attained below r reac . 

We denote by /? reac the heliocentric distance where T - T leac , 
which is given by: 



R„ 



I 3KmM 2 (GM ) i/2 

\1287r 2 cr SB ^BQ'r i 5 eac 



2/9 



= 1 AUM' 



4/9 



io- 



3\2/9 



0.5 m 2 /kg a 



500 K 



10/9 



(7) 



with M_g = M/(10 8 M Q • a : ), in the viscous-dominated 
regime. The mass accretion rate where R reac enters the 



irradiation-dominated regime (see equation (C.4i) is in our 
model: 



l287T 2 o- SB k B a(f T T ) 
3A-mQ^r r 4 eac 



9 ^ 



1/2 



/ AC 2/1 \ 1/2 / c \9/2 

4xl<r»M B .a- 1 ' " °- 5m /kg ^ lfT 



10" 



0.5 



500 K 



(8) 



which is a very low value (compared e.g. to mass accretion 



rates reported by Hartmann et al. 1998 1 so the disk is likely 
to have largely dissipated by that time. Moreover, it may be 
seen from equation (C.4i that then R leac < 0.1 AU, most likely 



inside the inner edge of the disk. So certainly at that time the 
validity of our disk model has broken down. We will find it 
convenient to take the corresponding time f leac as the final time 
in our calculation. 

Outside R Kac , the surface density of equilibrated water 
£H 2 o,eq i s governed by the same equation as total water, i.e. 
equation (|6jl with 2 H2 o replaced by S H ,o,eq- A stationary solu- 
tion to this equation, expected to be attained by f v ; s , is given by 
(Clarke and Pringle 1988; S tevenson|1990 see also Appendix 





\ R ) 



3Sc R /2 



(9) 



D/H of nebular water in a steady disk 
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Figure 2: D/H ratio (expressed in parts per million (ppm)) of water as a function 
of heliocentric distance in a steady disk, for three values of the radial Schmidt 
number Scr. In the "reaction zone", delimited by a vertical dashed line ("equili- 
bration line"), equilibrium isotopic fractionation with hydrogen gas is assumed 
; beyond, mixing with cometary water controls the D/H ratio (equation {TTJ). 
We have taken M = 1(T 8 Mo/a, a = 1(T 3 , K = 0.5 m 2 /kg, fy = 0.5. The 
smaller the Sen, the more efficient radial mixing is. The position of the snow 
line (where water condenses) is indicated by a vertical dotted line: water is 
gaseous inside the snow line and solid outside. 



where we have introduced the radial Schmidt number 

a 

Sc R = —. 
or 

We will henceforth adopt this solution. In other words, we as- 
sume that the radial distribution of equilibrated water under- 
goes a quasi-static evolution as the mass accretion rate (and thus 
R leac ) decreases. 

In this case, the D/H ratio of the total nebular water is given 

v3Sc 8 /2 



(10) 



(R,t) 



R 



eacWr 



(11) 



lh \\ 1L lh 

Thus, D/H is a monotonically increasing function both of he- 
liocentric distance and time (as over time M and hence 7? reac 
decrease). This is plotted in Fig. |2]and[3] 

2.4. Accretion and delivery to Earth of water-bearing chon- 
drites 

At this point, we have wholly prescribed the isotopic and 
transport properties of nebular water in our model. We have 
yet to relate this nebular water to the D/H distribution of chon- 
drites. 

We first need to make a prescription for coagula- 
tion/shattering, which determines the rate at which water is in- 



corporated in chondritic bodies. Similarly to 



(see also Ciesla and Cuzzi 2006 Weidenschil 



Cassen 



ing 



2004), we 



posit a collision term of the form 

£h 2 o 



6{R - tfcond(f)), 



(12) 



3 We use the fact that D/H<sc 1. 



D/H of nebular water in a steady disk 
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Figure 3: Same as Fig. [3]but for a fixed Sc«=0.3 and with varying mass accre- 
tion rates. 



with the coagulation timescale t C0iis (R) taken to scale like the 
local orbital timescale (Or 1 ), 6 the Heaviside function defined 
by 



0(z) 



Oifz < 
1 if z > 



(13) 



We thus consider that water is accreted solely as ice. 

We finally need to prescribe the delivery of accreted material 
to the Earth as a function of the heliocentric distance of initial 
agglomeration. While this step encompasses a variety of pro- 
cesses such as drift of meter-sized boulders, orbital evolution 
of parent bodies and eventually ejected meteoroids, we will be 
content in assuming a uniform probability of delivery of mate- 
rial to Earth up to a maximum initial heliocentric distance R mdx . 
In other words, the distribution we are calculating will be rep- 
resentative of water ice accreted inside R max . If the asteroid 
main belt accreted in situ, R mlLX could be taken to correspond 
to its outer edge, near 3 AU. If, on the other hand, significant 
redistribution of planetesimals occurred, e.g. during the "Grand 
Tack" studied by Walsh et al. ( 201 ^ max could be consider- 
ably larger. 
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3. Results 

Under the above hypotheses, meteoritic water, having been 
accreted from a range of heliocentric distances (between the 
snow line and R mlix ) and at various times, will necessarily ex- 
hibit a range of D/H ratios and one can define a probability dis- 
tribution function (PDF) for that quantity. The derivation and 



the expression of the PDF are presented in Appendix C and 
plots of it are shown in Fig. [4] It is notable that they are inde- 
pendent of a, to, ? C oag(l AU) 6h 2 o an d K an d depend on fj, T coa &, 
Treac, (D/H)/„ (D/H)/, R mSlX and Sc^?, of which only the latter two 
are considered free parameters. 

Distributions exhibit a minimum cutoff value (D/H) m j n q (see 



equation (C.8 1). This minimum value of the distribution stems 
from the constraint that the accreted water be condensed, and 



Figure 4: Probability distribution function of D/H calculated for different values 
of the radial Schmidt number Sen, assuming a maximum heliocentric distance 
of accretion R max of 3 AU (Top) and 10 AU (Bottom). The lower the Scr (that 
is, the larger the diffusivity), the more the PDF is peaked toward low values, 
while the larger the S max . the more the PDF is skewed toward heavy isotopic 
compositions. 
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Figure 5: Probability distribution function of D/H values in bulk carbonaceous 
chondrites (CI, CO, CV, CM) compared to the theoretical PDF adjusted to fit 
the minimum and maximum of the distribution. The comparison is warranted 
only if isotopic exchange with organic matter on the parent body is negligible. 
Data from|Kemdge|(l985);|Kolodny et a;.|p980);|Boato|(T954);|Robert and 
Epst eiii|jl982);|McNaughton et al.\\\9?,2y , Yang and Epstein|(1983) ; |Pearson 
et fl/.|p001) ; |Alexander et fl/.H2012bT . 



corresponds to the hydrogen isotopic composition of water at 
the snow line. It is thus governed by isotope exchange kinetics 
(which determines the composition and origin of equilibrated 
water) and radial transport to the snow line (which determines 
the proportion of equilibrated water there) alike. It is notewor- 
thy that the D/H ratio at the snow line does not evolve with time 
(so long it is in the viscous-heating dominated temperature re- 
gion), as it only depends on the ratio T con< i/T le - dc (see equation 

As is apparent on each panel of Fig. |4j low values of 
(D/H) m ; n o are associated with efficient radial mixing (enabling 
a high proportion of equilibrated water at the snow line), that is, 
a small radial Schmidt number, through the relationship (using 
equation (|C.8|l): 



Sc R 



„ ] n ( (D/H) A -(D/H), 
3 \ (D/HJ/j-fD/H),™ 

5 



ln(H 

V 1 cond / 



(14) 



For our fiducial parameters, if we adopt (D/H) m ; n o = 120 x 1CT 6 
from the observed PDF of D/H in carbonaceous chondrites, we 
obtain Sc s =0.2. 

A low value of Sc# is also needed if one is to recover 
the positive skewness of the observed PDF (with a negative 
slope over most of the range of D/H values). Indeed, equa- 
tion \C3\ indicates that the PDF has an overall dependence 
in [(D/H)/, - (D/H)] 1/(3Scs)_I , notwithstanding "second-order" 
modifications imposed by the constraint R con d < R < R mdx and 
the changes in the temperature regime. This imposes an over- 
all decreasing trend if Scr < 1/3^] Efficient transport indeed 
means that the D/H ratio will be close to that at the snow line 



4 Note that this inequality does not depend on the prescription we have 
adopted for M(t) but does depend on the power law for the irradiation- 



for a significant fraction of the outer disk, hence the dominance 
of this value in the PDF. 

By comparing the two panels of Fig. |4] it is also apparent that 
the PDF is also more skewed toward higher D/H ratio if R mdx 
increases, as is expected from the monotonic increase of D/H 
with increasing heliocentric distance in our model. The maxi- 
mum value (D/H) max j (see equation ( C. 1 1 1) is related to R mdx 



through the relationship (injecting equation ( 14 1 into equation 
dCTT) ): 



R n 



frTo 



exp 



f 10 ln (fc) ln ( 



(D/H),,-(D/H), 
(D/H^-fD/HW 



01 



ln( 



(D/H),, -(D/H), 
(D/H),,-(D/H)„ 



iii.O / 



(15) 



If we adopt (D/H) n 



= 230 x 10 , one obtains R mdx = 6 



AU. R mdx as expressed above is however quite sensitive to the 
assumed parameters, and as such our calculation would not con- 
clusively discriminate between in situ formation of the asteroid 
main belt and widespread redistribution (Walsh et al. 2 01 
Note that the evaluations of both Sc# and ^ max in equations 
( 14 1 and ( [15| are independent of our prescriptions of accretion 
of solids or the functional form of M(t). 

In Fig. [5] we have plotted the theoretical PDF with the Sc R 
and R mdx evaluated above (which by construction adjust the 
PDF to the minimum and maximum observed values and the 
histogram of the latter (observed) data). The shape of the ob- 
served PDF is reproduced qualitatively (though not quantita- 
tively), with (i) a peak near the lower end of the distribution 
and (ii) a tail toward high D/H ratios. The peak corresponds 
to D/H ratios near the snow line, which dominate the distribu- 
tion because accretion is most efficient there (because of shorter 
dynamical timescales and higher surface densities) than further 
from the Sun and also because efficient outward transport has 
almost homogenized the D/H ratios to the snow line value over 
an extensive region beyond it. The tail corresponds to water ice 
accreted some distance beyond the snow line, until the maxi- 
mum heliocentric radius sampled. 

However, the match is not quantitative, as the observed PDF 
is more peaked (around D/H = 150 x 10~ 6 ) than the theoretical 
prediction. Certainly, the simplifications used in the model (e.g. 
constant a, prescription of delivery to Earth etc.) prevent it from 
yielding realistic PDFs and only allow a proof-of-concept use. 
We shall now discuss implications of these results as well as 
assess the limitations of our treatment. 

4. Discussion 

4.1. Implications 

From the above results, it appears that the distribution of D/H 
ratios of carbonaceous chondrites can be accounted for in a sce- 
nario of isotopic exchange with hydrogen gas, with cometary 
D/H values prevailing at large heliocentric distances, under the 
condition that the radial Schmidt number Sc« be small (< 0.3). 



dominated temperature regime and that of the coagulation timescale: for 
T oc R-i and r coag oc R a , the 
(going back to equation jC.lj ) 
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The peak value of the carbonaceous chondrite population would 
be essentially dictated by the isotopic composition of water near 
the snow line, where accretion of condensed water would be 
most efficient. 

The low values inferred for Sc« would be consistent with hy- 
drodynamical turbulence (e.g. Pnrin[l990 l Dubrulle and Frisch 
1991 [ |Gail||2001|> — e.g. the 0.176 value measured by |Lath~ 



rop et al. ( 1992 1, although more experimental (and numerical) 
data would be desirable. It would not be consistent, however, 



with magnetohydrodynamic (MHD) turbulence ( Johansen et al. 



2006). MHD turbulence, mainly driven by the magnetorota- 



tional instability (MRI; Balbus and Hawley 1998 1, is widely 



believed to be the driver of transport in accretion disks. How- 
ever, there must be a region in the disk, referred to as the dead 
zone, which is too cold and dense for the MRI to operate ( Gam- 
|mie|1996) l. In the dead zone, turbulence should be of low inten- 
sity and hydrodynamical in nature. Thus our results are consis- 
tent with the idea that the planetesimals sampled by chondrites 
formed in this purported dead zone. In fact, ionization fraction 
calculations have shown that the dead zone would encompass 
the whole planet-forming region (from ~0. 1 to a few 10s of AU, 
e.g. |Bai and Goodman 2009 ) and indeed a dead zone would be 



a most favorable environment for planet formation (Terquem 
2008). Moreover the existence of a dead zone has been found 
to enable the preservation of chondrite components for a few 
Ma in the disk ( Jacquet et q/.|20lT) . 

The low value of Sc« would not only have enhanced out- 
ward transport of equilibrated water, but also that of higher- 



temperature components as well (see also Bockelee-Morvan 



et a/.|2002] >. This could account for the high-temperature min- 



erals, including calcium-aluminum-rich inclusions (CAIs) and 
chondrule fragments identified in comet Wild 2 samples ( |Zolen-| 
sky et al. |2006 Bridges et a/.|2012"] l and also the relatively high 



abundances of CAIs in carbonaceous chondrites (see |Jacquet 
\et a/.|2012| ) 

If mixing was as efficient as inferred above, the D/H ratio 
could have been significantly lower than (D/H)/, even at a few 
tens of AUs, depending on time (see Fig. G}. Thus, comets 
may have formed with D/H ratios lower than (D/H)/, (the limit 
at large heliocentric distances), which would account for the 



emerging range in measured D/H ratios of cometary water ( |Har- 



togh et al. 201 1 Bockelee-Morvan et al.\20V2\, consistent with 



the idea of an astero id-comet continuum ( |Gounelle et al. |2 008 ; 



Briani et al. 



2011 



^7 1 ' 

Even the 300 x 10 value adopted for 
(D/H)/, on the basis of measurements of known comets could 
then actually be a lower bound of its real value. 

4.2. Modeling caveats 

Of course the validity of the conclusions of this study can 
only be as good as that of its underlying assumptions. 

First, the quasi-static approximation we have used requires, 
as we recall from Section |2~Tj that the viscous timescale t v i s (R) 
be short relative to the time of chondrite accretion and that of 
evolution of the mass accretion rate. In a dead zone, whose exis- 
tence we have inferred above, a would likely be small, around 
10~ 4 (e.g. |Fleming and Stone 



2003 



Ilgner and Nelson 



2008 



D/H of nebular water in a steady disk 




10 20 30 

Heliocentric distance in AU 



Figure 6: The D/H ratio as a function of heliocentric distance of accretion, 
for the radial Schmidt number inferred from this study (0.2) at a fiducial time 
where the snow line lies at 3 AU from the Sun. We have taken k = 0.5 m 2 /kg 
and a = 10~ 4 (as appropriate in a dead zone). To realistically extrapolate to 
the comet-forming regions, we have here taken into account the effect of radia l 
drift by considering millimeter-sized grains [p s a = 1 kg/m 2 ; see | Appendix A| . 
The D/H curve without this effect is also shown for comparison (dotted line). 



be longer than our nominal evaluation in equation ([3]), but at 
R = 3 AU, it would be ~ 1 Ma, still shorter than the accretion 
time of chondrites after the start of the solar system (~l-5 Ma, 
see e.g. |Villeneuve et al. ( 2009| l; Connelly et al. ( 2012) )). Ad- 
mittedly, the timescale constraint would be quite marginally sat- 
isfied, so that, conceivably, future re-examination of this prob- 
lem in time-dependent simulations of disks with dead zones, 
with self-consistently evaluated a (see e.g. |Zhu et aT\\10\0\ ), 
could unveil new interesting effects. Nonetheless, this would 
not affect our conclusion that MRI-turbulent disk models, with 
higher a — hence largely satisfying the timescale constraint — 
and Scr, cannot reproduce the D/H ratio of chondrites and our 
inference of the presence of dead zone thus appears robust in 
this respect. 

Our steady-state disk model neglects photoevaporation dur- 
ing the bulk of the disk's lifetime. However, if, as argued by 
Desch ( 2007 | l, the disk was subject to intense photoevapora- 
tion from the outset, the net gas flow could have been outward, 
and equilibrated water would have reached the snow line al- 
most undiluted; i.e. water at the snow line would have a D/H 
ratio close to (D/H)/ as 40 x 10~ 6 . If such isotopically light bulk 
preaccretionary water compositions were to be found in the fu- 
ture, this would be a possibility worth investigating ; however, 
to date, such isotopic compositions have not been measured. 

In this work, nebular water has been assumed to be dynami- 
cally coupled to the gas. Depending on the fragmentation and 
the bouncing barriers (e.g. |Zsom et~al. 2010 Birnstiel et al. 



2012| l, it is however conceivable that ice-bearing particles grow 
to millimeter-size and beyond so that drift due to gas drag 



Oishi and Mac Low 2009 1 |Turner et al. |20 1 1, so that f v i S would 



would be significant (see equation ( A. 5 1) before they are incor- 
porated in planetesimals. This would likely lead to significant 
enhancements of water abundance inside and at the snow line, 
and depletions further outwards (e.g. Cuzzi and Zahnle 2004 
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Ciesla and Cuzzi 2006 ). From equation (A.81 derived in Ap 



pendix A it can be seen that the fraction of equilibrated wa- 
ter would be lower, and thus D/H higher, than under our tight 
coupling calculation, because of increased import of "cometary 
water". One example of this effect is shown in Fig. [6] us- 
ing the parameters inferred from this study. While the effect 
might be limited for the carbonaceous chondrites which we 



have focused on — for which Jacquet et al. (2012 1 inferred that 
millimeter-sized components were not significantly drifting — , 
it would certainly accelerate the convergence of the D/H ratio 
toward its asymptotic value in the comet-forming region (for 
Sr > 1, see Appendix Al. This could account for the apparent 



carbonaceous chondrite/comet hydrogen isotopic dichotomy, in 
the sense that relatively little material with intermediate com- 
position would then exist (but does exist, as emphasized at the 
end of Section |4. \\ . However, regardless of the actual quanti- 
tative importance of radial drift, it must be noted that inasmuch 
as nebular water inside the snow line would be in vapor state 
— even meter-sized solids would lose their water by vaporiza- 



tion before they can travel much inward of it (Cuzzi and Zahnle 
2004 ) — and thus continue to be tightly coupled to the gas, the 
transport of equilibrated water until the snow line would be un- 
changed (see equation A. 8 1 — save for some transition period 
preceding establishment of the quasistatic regime. Therefore, 
our estimate of Sc« from the minimum D/H of the distribution 
(established at the snow line) would be unaffected and is thus 
robust in this regard too. 

We have also argued that accretion was inefficient, allowing 
us to ignore feedback on the transport of nebular water. If this 
approximation were not to hold, rapid accretion of planetesi- 
mals at the snow line could progressively deplete the inner so- 
lar system in water (e.g. |Ciesla and Cuzzi||2006| l. Although, 
again, this would likely not affect the D/H ratio at the snow line, 
the effects on the shape of the PDF have yet to be determined 
by dedicated numerical simulations including both coagulation 
and shattering (see Yang et al.\l0\2) . 



4.3. On the interpretation of D/H chondrite data 

In this work, we have compared our PDF of the isotopic com- 
position of water to bulk chondrite compositions (see Fig. [5j, 
and, in doing so, implicitly assumed that the latter reflected the 
composition of pre-accretionary water now locked in hydrated 
silicates. However, hydrated silicates are not the only contrib- 
utors to the hydrogen budget of carbonaceous chondrites, as 
the latter also contain organic matter. The rationale for making 
this identification nonetheless is that inasmuch organic matter 
does not account for more than 10 % of the hydrogen, at their 
presently measured D/H ratios (mostly < 400 x 10~ 6 — still ex- 
cluding CR chondrites — ; Alexander et al. (2010 1), the incurred 
shifts to the isotopic composition would be small (< 20 x 10 -6 



Robert (2006 )) compared to the observed range of D/H ratios. 



However, this ignores the possibility of isotopic exchange be- 
tween organic matter and water on the parent body, in which 
case the organic matter could have been much more D-rich ini- 
tially than presently measured. In fact, | Alexander et al. ( 2012b] ) 
showed that the D/H ratio of bulk carbonaceous chondrites cor- 
related positively with the C/H ratio, and interpreting this as a 



mixing trend, extrapolated an initial D/H ratio for organic mat- 
ter comparable to that of CR chondrites (« 700x 10~ 6 ), and con- 
sequently, a lighter isotopic composition for pre-accretionary 
water, e.g. (86.5 ± 3.5) x 10" 6 for CM chondrite^] In that case, 
while (direct) comparison of the theoretical PDF with the his- 
togram plotted in Fig. [5]would no longer make sense, (D/H)^,, q 
would be constrained to be lower than this latter value, and our 
formalism would thus constrain Scr to be < 0.1 (see equation 



14 1), that is, a yet more efficient radial diffusion would be in- 



dicated. 

It must be cautioned, though, that in situ observations of CI, 
CM and CR chondrite matrices have revealed little evidence 
of isotopic exchange between the isotopically heterogeneous 



organic matter and the intermixed hydrous minerals ( Remusat 



\et aZ.|[2010] l and organic matter in the least aqueously altered 
CM chondrite Paris has an hydrogen isotopic signature indis- 
tinguishable of that of the other CM (Remusat et al. 201 1 1. One 



possible explanation for the D/H - C/H correlation of|Alexander 



et al. ( 2012b I could be that the cometary water endmember was 



somehow coupled to D-rich organic matter in the disk, consis- 
tent with the high carbon contents of cometary grains (e.g. |Woz^1 



niakiewicz et a/.|2012| ). In that case, the mixing trend would not 



actually extrapolate to the composition of the organic matter 
but to that of the resulting composite (water+organics) C- and 
D-rich endmember. In this picture, the hydrogen isotopic signa- 
tures of organic matter and hydrated silicates could be actually 
largely pristine, as inferred by Remusat et al. (2010). What- 
ever that may be, and whether D/H ratios of bulk carbonaceous 
chondrites truly reflect the D/H ratios of preaccretionary water 
or only provide upper limits, we stress that low Schmidt num- 
bers are robustly required in the framework of our formalism. 

While the present work has mainly focused on carbonaceous 
chondrites, it is noteworthy that non-carbonaceous chondrites 
are generally richer in deuter ium than the former (e.g. |Robert| 



2006 



Alexander et aZ.||201 2bp1 As the D/H ratio of nebular 



water increases with time in our scenario, this could imply that 
they accreted later than the former, as suggested by |Jacquet| 
\et al.\ ( |2012| ) based on the modeled redistribution of chondrite 
components in the disk. Conceivably, radial drift, as alluded to 
in the preceding subsection, which would be most pronounced 
as the coupling of the grains with the less dense gas would 
be looser then, might have contributed to significant enhance- 
ments of the D/H ratios for them. In fact, Jacquet et al. ( 2012| ) 
specifically inferred that millimeter-size components were sig- 
nificantly drifting for these chondrites, and from equation ( |A.8[ >, 
the nebular water should be isotopically close to cometary val- 
ues, which appears to be the case for clays from unequilibrated 
ordinary chondrites (e.g. |Robert|2006[|Alexander et a/.|20 12b ). 



5 This would incidentally explain some low-D/H point measurements in LL3 
cho ndrites IDeloule e t a/. 1 199 8). 

^Alexander et al.\ |2ul2a} suggest that the high D/H ratios of these mete- 
orites could be due to oxidation of metallic iron by water producing isotopi- 
cally light hydrogen gas escaping from the system. However, while this may 
have affected metamorphosed chondrites I Alexander et al. 20 1 McCanta etal. 
2008) , the most unequilibrated ordinary chondrites also show these high values, 
with considerable heterogeneity consistent with the pristinity of this signature 
jPiani|2012} . 



A later formation of the non-carbonaceous chondrites would 
avoid the alternative possibility of an accretion further from the 
Sun than carbonaceous chondrites, which would run counter to 



the observed distribution of asteroid classes (e.g. Burbine et al. 
2008). A later formation time could also account for the high 
D/H values of CR chondrites too — consistent with their young 
chondrule Al-Mg ages (Kita and Ushikubo 20121 — , which 
could mark an intermediate status between the other carbona- 
ceous and the non-carbonaceous chondrites — along with their 
rather low refractory inclusions abundance, their limited re- 
fractory element enrichment relative to Mg, and their relatively 
heavy oxygen isotopic composition ( Scott and Krot|2 003). 



5. Summary 

We have considered a simplified analytic model for water 
transport and isotopic evolution in an evolving disk based on 
the following main assumptions: 

(i) The regions of interest of the disk are approximated by a 
stationary model, with a constant turbulence parameter a. 

(ii) Accretion is inefficient and does not significantly affect the 
transport of nebular water. 

(iii) Water far from the Sun is assumed to have the same D/H 
ratio as cometary water, and any batch of water having 
experienced temperatures above a "reaction temperature" 
r reac has had its D/H reset to a fixed value (D/H)/. 

(iv) Accretion is modeled by a coagulation timescale scaling 
with the orbital period and water is incorporated as ice 
only. A uniform probability of delivery to Earth is applied 
to agglomeration locations up to a maximum heliocentric 
distance R max . 

In this model, the D/H ratio of water is an increasing function 
of time and heliocentric distance. After integration over time, 
the probability distribution function (PDF) of the D/H ratio in 
accreted water is found to depend essentially on R max and (most 
sensitively) on the radial Schmidt number Sc s = a/S^ with <S S 
parameterizing the diffusivity. The minimum cutoff of the PDF 
is determined by isotopic composition of water at the snow line, 
while the isotopically heavy tail is dictated by R mdx . 

It appears that the model is able to broadly account for 
the observed PDF in carbonaceous chondrites if low values of 
Scr (around 0.1-0.3) — i.e. efficient outward diffusion — are as- 
sumed in order to reproduce the low D/H values of most chon- 
drites and the positive skewness of the observed distribution. 
This would be most consistent with hydrodynamical turbulence 
as expected to prevail in the dead zone of the protoplanetary 
disk. Efficient outward diffusion would also have enabled the 
transport of high-temperature minerals to comets. The high 
D/H ratios in CR chondrites and non-carbonaceous chondrites 
could indicate an accretion later than most carbonaceous chon- 
drites. 

The effects of radial drift and higher accretion efficiencies 
on the transport of water and its hydrogen isotopic composition 
have yet to be investigated. 
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Appendix A. Diffusion of equilibrated water in disks 

In this appendix, we calculate the steady-state profile of neb- 
ular water concentration and its equilibrated fraction. To allow 
discussion of the tight coupling assumption in Section [4] we 
take into account the finite size of ice-bearing particles beyond 
the snow line, and thence the effects of gas drag, so that their 
velocity becomes, ignoring any feedback of the solids on the 
gas (e.g. Birnstiel et a/.|2010 i: 



1 



vr 



a 



i + St 2 \ s/? 1/2 OR 



rdP 
pdR 



(A.l) 



with v = ac 2 s I Q, the turbulent viscosity, r the gas drag stopping 
time, St = Qr the Stokes number, and P and p the gas pressure 
and density, respectively. 

In steady state, the mass accretion rate of nebular water 



M H2 o = 2nR (-S^oV* + £>«X^ f^jj , 



(A.2) 



is constant. We have introduced the diffusion coefficient modi- 
fied by finite particle size as ( |Youdin and Lithwick|2007] l 



D R = 



1 + St 2 ' 



(A.3) 



Equation \A.2) may be viewed as a first-order ordinary differen- 
tial equation in Z H ,o/£- It is noteworthy that the corresponding 
homogeneous equation is the equation governing the transport 
of equilibrated water, since its concentration vanishes at infin- 
ity, so that: 



exp 



J d r 



exp(f^dK') 



(LvR 1 / 2 ) 



3Sc„ 



(A.4) 



where the second proportionality relationship assumes that Sc« 
is radially constant and we have coined Sr = St/6R which is 
a measure of gas-grain decoupling (Jacq ueTef <//.|[2~()l2). In 
the Epstein drag regime, for spherical grains of density p s and 
radius a (averaged over the size distribution), we have: 



S R = - 



n p s a 



2 2<5 S 



37t 2 Sc s p s ac- s 
2 MQ. 



0.1- 



p s a 



M_ % \lkg/m : 



)(— ),(A.5) 
/\300K/' V ; 



where the last two equations pertain to a steady disk as assumed 
in the main text. Note that the normalizing value p s a = 1 kg/m 2 
corresponds to millimeter-sized grains (the typical size of non- 
matrix chondrite components). 
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By requiring that £ Hn o/£ does not diverge at the disk's inner 
edge (whose heliocentric distance we denote by R m , taken to be 



zero in the main text), equation ( A. 2 1 may be then integrated as: 

^Uexpff *dtf) f exp(- 
2 Ik^* /J« m P \ J«, n D R ]2nR'XD R 

For 5s <K 1 (in particular inside the snow line where S r — 0), 
this is a constant , and for Sr » 1 » St, it falls off as 1/Sr. 
Then, the fraction of equilibrated water is: 



(A.6) 



-H 2 0,eq 

£h 2 o 



fexpf-r^")^- 



(A.7) 



Enforcing that the left-hand-side be unity at R — R le:ic , this be- 
comes: 



-H 2 0,eq 

£h 2 o 



3Sc R 



-) 



(v^-v^) 3sc " r (v^- 

( C R ' „ dlnP A, „ 9 x d/?' 



For Ss « 1, this amounts to equation (|9]l for R[ n <k /? reac < 

Appendix B. Kinetics of water-hydrogen isotopic exchange 

Ignoring transport, and given that D/H <c 1 and H2O/H2 <s 
1 , the equation governing the evolution of D/H of water may be 
written as ( |Lecluse and Robert|1994| i: 

(T 1 )- I 77 1 1 (B.D 




3Sc«-l 



d /D 
df\H 



r(T)«H 2 



H 



with the number density of hydrogen molecules, (D/H) eq 
the equilibrium values (see e.g. Richet et al.\\911\ and the rate 
constant 



k~(T) = 2 x 10- 28 expl-^^ )m 3 /s. 



(B.2) 



Given that, in our steady-state model, the midplane number 
density can be expressed as a function of temperature (assumed 
to be vertically constant) as 

1/3 / 27^1 



V2tto 



16/ 3 
9 \Ma 



nk„m , 



\ K I 



the characteristic equilibration timescale resulting from equa- 
tion (IB. lb is: 



1 



& (r)«H 2 

exp (™) / 



= 0.1 Ma 



,11/6 



I 10 ~ 



1/3 



0.5 m 2 /kg 



2/3 

(B.4) 



This is a sharply decreasing function of temperature. While, 
at high temperature, f eq is short compared to the transport 
timescale f v ; s so that water vapor and hydrogen gas can be con- 
sidered in equilibrium, at lower temperature (at larger heliocen- 
tric distances), f eq becomes long compared to f v is and isotopic 



Figure C.7: A schematic diagram to help visualize the calculation and the dif- 
ferent quantities involved. The abscissa is the D/H ratio (x) and the ordinate is 
the time f, or, correspondingly, the mass accretion rate M. The snowflakes mark 
(rA(8>)nain where water is condensed and those colored in deep blue mark water 
ice inside heliocentric distance R max , which contributes to the PDF. Their area is 
bounded by the thick black curves (x mm and .r max ), while the dashed line marks 
the boundary between the viscous-heating and irradiation-dominated regimes 
(for the location corresponding having the D/H value of x at time t). With the 
passage of time, the average D/H ratio of water increases. The blue line rep- 
resents the time interval of integration on the right-hand-side of equation jC.3) 
for a given value of x. 



exchange is essentially quenched. The temperature r reac mark- 
ing the transition between the two regimes can be determined 
by setting: 



^eq(-^reac) — ^vis(-^reac)- 



(B.5) 



Water originating from inside R Kac ("equilibrated water") will 
then essentially have the isotopic composition it had when it 
last equilibrated with hydrogen gas, i.e. at T — T reRC , that is, its 
D/H ratio will be (D/H)/ =(D/H) eq (r reac ). 



Numerically, equation (B.5 1 can be expressed as 



exp 



(51TOK\ 

V reac / 



5 x 10 3 no 



r 5/18 
reac,K 



M 



1/9 



3 v 13/9 



0.5 m 2 /kg 



,5/9 



(B.6) 



The left-hand-side being a sharply decreasing function of T le!LQ , 
the sensitivity on M is completely negligible and that on a is 
also fairly weak. We shall thus adopt the solution of this equa- 
tion for a = 10~ 3 , M = 10~ 8 M /a, k = 0.5 m 2 /kg, which is 
TVeac = 500 K. Adopting a protos olar D/H ra tio (for the H2 gas) 
of (20±3.5)xl 0~ 6 (|Geiss and Gloeckler|2003t , the fractionation 
factor given by|Richet et al.\\\91l\ yields (D/H)/ = 40 x 10~ 6 . 



Appendix C. Calculation of the PDF of D/H 

In this appendix, we detail the calculation of the PDF of the 
D/H ratio of meteoritic water. We will here denote D/H by x for 
the sake of legibility of the equations. 
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With the model set up in Section [2] the mass of meteoritic 
water with D/H ratio between x and x + dx is 



(* 1[cac dR 

S coU (R(x,t),t)2nR(x,t) — dx9{R n 
Jo ox | f 



■ R(x, t))dt 
oc/(x)dx(C.l) 



where f(x) is the normalized probability density function of 
bulk chondrites in terms of D/H ratio, t — is taken to corre- 
spond to the time where R con d = R m ax, i-e. the earliest possible 
time where water can condense and be accreted inside R mm , 
corresponding to the mass accretion rate 



M = 



( 1287rV SB ^r 5 on / 
3/<mfi(fl max ) 3 



1/2 



= 1.5x 10" 7 Mo.a _1 



(-^cond \ / ^max \ 
170 k) Iioau/ 



a 0.5 m 2 /kg 



io- 3 

9/4 



1/2 



(C2) 



Using the assumed constancies of e^o, D.t coag (R), a, Sc« and 
using equation we have: 



fix) 



Rte 



M(t) 



(x h - 



oit)T(R(x, t), t) 
0(R mm - Rix, t))6(R(x, t) - R cond (t))dt 



,1 



••/.\'/(/)\ 

"So"/ 



7/9 



dr 



(C3) 



with t m i n (x) and ? ma x( x ) determined by the conditions R < R max 
and R > R con d, respectively, M\, r (x) the mass accretion rate 
for which x(R bl ) = x, where R br is defined as the heliocen- 
tric distance of the transition between the viscous-heating- and 
the irradiation-dominated regimes and is given by (see equation 



^br. 



M 2 D 3 1 



mn 2 a- SB k B a (/ r r ) 5 J 



1/2 



M_ 



10 



-3 \ 1/2 



(C4) 



/5/ 2 \0.5 m 2 /kg a 

with f2o = Q(l AU). The "min(...)" factor in equation ( |C.3| l is 
unity when R(x, t) is in the irradiation-dominated regime. 

Before proceeding to the result, we coin a few additional no- 
tations and then proceed to the integration. Fig. C.7 may help 
the reader to visualize the situation in the D/H - time space. 

For a given mass accretion rate, the minimum and maximum 
value of D/H dictated by the conditions R con d < R < R m <tx are: 



XmmiM) = 



5Sc«/3 



Xh 



x h - (x h - *;) (§>ff) if T bl < r cond 
_ (XI _ x ,)((^) 3 ( 3 *"^ ) lf3 ^' 

(X h Xl)\\f TT J \u^cr SB k Ba T^) 



and 



x max(M) - X h - (x h - X[) 



3Sc s /2 



1 



(C.6) 



The D/H ratio at R - R br (for M > M leac ) is: 

fl2,S^ ( r SB A, () (/ / '/n)"^ S> 



Xbr(M) = x h - (x h - Xi) 



3 K mM 2 Q 3 Ti. dc 



(C.7) 



To all these, a "0" subscript will be added if evaluated for 
M = Mo and "f ' for M = M rea c- Hence, 



-^min,0 — -^max,0 — %h CVft X/) 
^irr(^max) 



- cond 



5Sc fi /3 



*br,0 = Xh~ (Xh - Xi) 



T K . 

9\5Sc ff /12 



T T 4 T^5 

V 1 reac 1 con( j 7 



= x A - (xh - xi) 



Tcond \ 

7^ 



3Sc R 



-^irr(^max) 



3Sc ft 



(C.8) 



(C9) 



(CIO) 



(C.ll) 



— 1/2 

where ri rr (/? max ) = frToR ' AU is the irradiation temperature 
at heliocentric distance /? max . 

We will also need the D/H ratio when 7?br = R m ax> which is 
given by: 



x maxbr — x /i ( x /i Xj) 



T'in^^max) 



5Sc s /3 



(C.12) 



and the ratio M reac /Mo between the final and the starting mass 
accretion rate: 



M„ 



T'irr^max) 



T^4 T^5 

V J reac J cond / 



1/2 



(C.13) 



With all these notations, we can express the result of the in- 



tegration in equation (C.3 i as: 



f(x)=C(x h -x) lK3Sc « ) - 1 



T 13 / ls T 5/w \x h -x 

reac cond 



13/(30Sc«) 



Xh - Xhr,0 
Xfi X 

M, 



2/(15Sc«) 



Mo 



1/9 



Six) 



(C.14) 



lZST^CTSB^flOfTi 

if ^cond < ^br < T, 
\3Sc fl 



for Xmin.o <x< x maxbr , and 

f{x) = C(X A - jc)V(3Sc s )-1 

(C.5) 



Xfo X 
Xfi -^min,0 

M r 



l/(6Sc«) 



M 



1/9 ^ 

Six) 



(CAS) 



Xh ~ixh-xi){^) " if T br > r reac for x maxbr < x < x max>f . 
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Here, C is a normalization factor given by 



C = 



Sc^a - *0 1/(3Scs) 12 



, r,„(Rmax) 1/2 r!; 
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JkrWnax) 5 ' 6 



Tlu r (/I T 5/9 T 5/18 T ll/9 7-11/9 I 

i cond i rcac ^cond^reac ' reac y 



and S (x) is defined as 



S(x) 



l/(6Sc R ) 



if x < x, 
if x > x n 



min,f 
in,f 



(C.16) 



(C17) 
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